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Abstract 

We present a general scheme and scale independent parameterization of two-body non-leptonic B 
decay amplitudes which includes perturbative QCD corrections as well as final state interactions in 
a consistent way. This parameterization is based on the Next-to-Leading effective Hamiltonian for 
non-leptonic B decays and on Wick contractions in the matrix elements of the local operators. Using 
this parameterization, and making no dynamical assumption, we present a classification of two-body B 
decay channels in terms of the parameters entering in the decay amplitudes. This classification can be 
considered as the starting point for a model-independent analysis of non-leptonic B decays and of CP 
violating asymmetries. We also propose, on the basis of the large N expansion, a possible hierarchy 
among the different effective parameters. We discuss the strategy to extract the most important 
effective parameters from the experimental data. Finally, we establish a connection between our 
parameterization and the diagrammatic approach which is widely used in the literature. 
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1 Introduction 



A quantitative description of two-body non-leptonic B decays in the framework of the Standard Model 
remains as an important chahenge for theorists. Simultaneously these decays play a decisive role in the 
study of CP violation and the determination of the Cabibbo-Kobayashi-Maskawa (CKM) parameters 
at i?-factories and dedicated -B-physics experiments at hadron colliders. The studies of these decays 
should also provide some insight into the long distance non-perturbative structure of QCD. 

The basic theoretical framework for non-leptonic B decays is based on the Operator Product 
Expansion (OPE) and renormalization group methods which allow to write the amplitude for a decay 
of a given meson B=Bd, Bg, B^ into a final state F=7nr, Kir, DK, KK, . . . , generally as follows: 

AiB ^F) = {F\nes\B) = '^Y.^r''C,{^i){F\QMB). (1) 

Here Weflf is the effective weak Hamiltonian, with Qi denoting the relevant local operators which 
govern the decays in question. The CKM factors V^^km ^j^^ Wilson coefficients Ci{fi) describe 
the strength with which a given operator enters the Hamiltonian. In a more intuitive language, the 
operators Qiin) can be regarded as effective vertices and the coefficients Cj(/u) as the corresponding 
effective couplings. The latter can be calculated in renormalization-group improved perturbation 
theory and are known including Next-to-Leading order (NLO) QCD corrections Q^Q- The scale fi 
separates the contributions to A{B F) into short-distance contributions with energy scales higher 
than n contained in Ci{p) and long-distance contributions with energy scales lower than /i contained 
in the hadronic matrix elements {Qi{fj,)). The scale n is usually chosen to be 0{mb) but is otherwise 
arbitrary. The //-dependence of Ci(//) has to cancel the /i-dependence of {Qi{n)) so that the physical 
amplitude A{B — > F) is /U-independent. Similarly, the renormalization scheme dependence of Ci{n) 
cancels the one of {Qi{fi)). It should be stressed that these cancellations involve generally several 
terms in the expansion (||). 

The great challenge for theorists is a reliable calculation of the matrix elements {Qi{n)) in QCD. 
Unfortunately, due to the non-perturbative nature of the problem, the progress towards this goal has 
been very slow and it is fair to say that no reliable calculations of the hadronic matrix elements {Qi{f^)) 
in QCD exist at present. In the case of two-body non-leptonic B decays, the {QiifJ-)) cannot in general 
be calculated from first principles in lattice QCD, due to the Maiani- Testa no-go theorem [^, and 
even the feasibility of a model-dependent estimate of these matrix elements by numerical simulations 
still has to be verified |^]. 

In view of this situation a number of strategies for the calculation of A{B — > F) have been used 
in the literature. The most extensive analyses have been done in the factorization approach [^-|^ in 
which the hadronic matrix elements {Qi{n)) are replaced by the products of the matrix elements of 
weak currents. The latter can be expressed in terms of various meson decay constants and generally 
model-dependent form factors. Thus the decay amplitude in this framework is given simply by 

Ai,n = ^l^^™ai,2(/i)(Qi,2(/i))F, (2) 

with 

Ol,2(/^) = Ci,2(//) + ^C2,l(/i) (3) 

and being the number of colours. 
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Here {Qi{^j))F denote the factorized matrix elements of the current-current operators Qi^2 given 
expHcitly in Section The indices I and II distinguish between the so-called class I and class II decays. 
Since in this approach the matrix elements are scheme- and ^u-independent, the resulting 

amplitudes have these dependencies, which are clearly unphysical. 

Two ways of remedying these difficulties have been suggested, which led to the concept of gen- 



eralized factorization [10|-||lJ]. In the formulation due to Neubert and Stech |T^], the //-dependent 
parameters ai^2(A*) are replaced by the and scheme-independent effective parameters af'l- The 
latter depend formally on Cj(/t) and non-factorizable contributions to (Qi(/x)) which are supposed to 
cancel the ^- and renormalization scheme dependences of ai{^). In this framework there is no explicit 
calculation of non-factorizable contributions and af^ are treated as free parameters to be extracted 
from the data. With the assumption of universality and the neglect of a class of non-factorizable 
contributions represented by penguin diagrams and Final State Interactions (FSI), two-body decays 
in this approach are parameterized by the two real free parameters 

The generalized factorization presented in is similar in spirit but includes more dynamics 



than the formulation in |10|. Here the non-factorizable contributions to the matrix elements are 
calculated in a perturbative framework at the one-loop level. Subsequently these non-factorizable 
contributions are combined with the coefficients Ci [jj) to obtain effective and renormalization scheme 
independent coefficients Cf^ . The effective parameters are given in this formulation as follows: 

1 1 

with analogous expressions for af^ (i = 3 — 10) parameterizing penguin contributions. Here A'^'^*^ 
is treated as a phenomenological parameter which models the non-factorizable contributions to the 
hadronic matrix elements. In particular it has been suggested in |]l^]-|13| that the values for 



extracted from the data on two-body non-leptonic decays should teach us about the pattern of non- 
factorizable contributions. 

A critical analysis of these two approaches has been presented by us in ref. [^]. In particular 
we have pointed out that the effective coefficients Cf^ advocated in |11, 13] are gauge and infrared 



regulator dependent. This implies that the effective number of colours extracted in [11, ^] also carries 
these dependencies, and therefore it cannot have any physical meaning. Concerning the approach 



of Neubert and Stech [10|, we do not think that this approach is most suitable for the study of 
non-factorizable contributions to non-leptonic decays. In particular, in the present formulation the 
contributions from penguin operators, penguin diagram insertions and final state interactions are 
not included. This implies, for instance, that a large number of penguin-dominated decays cannot 
be properly described in the present formulation of this approach. Another important issue in the 



analyses performed in generalized factorization is the dependence on the form factors used. See ref. [16| 
for a detailed discussion of this point. 

Another, more general, approach to non-leptonic decays is the diagrammatic approach, in which 
the decay amplitudes are decomposed into various contributions corresponding to certain flavour- 
flow topologies which in the literature appear under the names of "trees" , "colour-suppressed trees" , 
"penguins", "annihilations" etc. |l^, Supplemented by isospin symmetry, the approximate SU{?>) 
flavour symmetry and various "plausible" dynamical assumptions the diagrammatic approach has 
been used extensively for non-leptonic B decays in the nineties. 

Recently the usefulness of the diagrammatic approach has been questioned with respect to the 
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effects of final state interactions. In particular various "plausible" diagrammatic arguments to neglect 
certain flavour-flow topologies may not hold in the presence of FSI, which mix up different classes of 
diagrams |1£]-[26|. Another criticism which one may add is the lack of an explicit relation of this 
approach to the basic framework for non-leptonic decays represented by the effective weak Hamiltonian 
and OPE in eq. ([l|). In particular, the diagrammatic approach is governed by Feynman drawings with 
W-, Z- and top-quark exchanges. Yet such Feynman diagrams with full propagators of heavy fields 
represent really the situation at very short distance scales 0{Mw^z,iT^t), whereas the true picture 
of a decaying meson with a mass 0{mb) is more properly described by effective point-like vertices 
represented by the local operators Qi. The effect of W, Z and top quark exchanges is then described 
by the values of the Wilson coefficients of these operators. The only explicit fundamental degrees of 
freedom in the effective theory are the quarks u, d, s, c, b, the gluons and the photon. 

In view of this situation it is desirable to develop another phenomenological approach based directly 
on the OPE which does not have the limitations of generalized factorization and allows a systematic 
description of non-factorizable contributions such as penguin contributions and final state interactions. 
Simultaneously one would like to have an approach that does not lose the intuition of the diagrammatic 
approach while avoiding the limitations of the latter. 



First steps in this direction have been made in refs. |2J, 26, In refs. [|24| , p?] ] some of the 
parameters of the diagrammatic approach have been, in the case of B ^ Kir, expressed in terms 
of matrix elements of local operators. On the other hand in ref. [^6| the amplitudes for B — > tttt, 
B Kir and B KK have been given in terms of diagrams representing Wick contractions of the 
operators of the effective Hamiltonian between the relevant hadronic states. 

Now the formulation of non-leptonic decays given in ref. |26] involves the explicit expansion of 
the decay amplitudes in terms of Wick contractions; therefore, with ten operators entering the basic 
formula (||) and several possible contractions, one ends up with over hundred different contributions, 
each of them being scale and renormalization scheme dependent. 

In the present paper we make the approach of ref. manifestly scale and scheme independent. 
In this context we introduce a set of effective scale and scheme independent parameters which are 
given as linear combinations of particular Wick contractions of the operators times the corresponding 
Wilson coefficients. This reformulation of the approach in [^] in terms of effective parameters results 
in more transparent formulae for the decay amplitudes than given in [26|, establishes some connection 
with the usual diagrammatic approach and is more suitable for approximations. Moreover we include 
additional Zweig-suppressed Wick contractions, not considered in the literature, and we study in 
addition to charmless final states also those including the charm flavour. 

The formulae for the effective parameters given here allow in principle their calculation in QCD by 
means of lattice techniques or other non-perturbative methods. Since this is not possible at present 
they have to be considered as free parameters to be determined from the data. 

It turns out that the full description of two-body S-decays requires the introduction of fourteen 
flavour-dependent effective parameters. With the help of large ideas and plausible dynamical 
assumptions one can argue that several of these parameters play only a minor role in two-body non- 
leptonic decays. 

The approach presented here allows for a general phenomenological description of non-leptonic 
decays which with more data could teach us about the role of non-factorizable contributions and 
about the flavour structure of non-leptonic decays. In order to be more predictive some symmetry 
relations, as S'C/(3) relations [28|, and dynamical input are needed. In this context we would like to 
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mention an interesting work on non-factorizable contributions to non-leptonic decays within the QCD 
sum rules approach ||2^. We will return to several of these issues in a subsequent publication. 

Our paper is organized as follows. In Section |2| we recall the complete effective weak Hamiltonian. 
In Section ^ we classify various topologies of Wick contractions. We identify fourteen flavour topologies 
of which only nine have been considered in the literature. The new topologies correspond to Zweig- 
suppressed transitions. They may play a role only in certain decays but strictly speaking they have 
to be included for consistency. 

In Section ^ we introduce the effective parameters in general terms. In Section ^ we derive by means 
of a diagrammatic technique the explicit expressions for these parameters including flavour dependence. 
In Section |6| we propose, using the expansion, a hierarchical structure for the effective parameters. 
In Section ^ we classify the and decays into suitable classes and we give explicit expressions for 
a large number of decays in terms of the most important effective parameters, following the hierarchy 
proposed in Section ^. In Section ^ we give analogous expressions for Bg decays. The contributions 
neglected in the analysis of Sections |^ and ^ are collected in Appendix |^ for completeness. In Section 
^ we discuss briefly strategies for the determination of some of the effective parameters from the data 



and in Section 10 we compare the present approach with the diagrammatic approach of refs. |17, 18 1. 
We end our paper with a brief summary and conclusions. A detailed application of this formalism to 
two-body B decays, in particular to CP asymmetries, will be presented elsewhere. 



2 Effective Hamiltonian 



The complete effective weak Hamiltonian for non-leptonic B decays is given by: 



n 



eff 



^ubVud 



+yubv:s 
+yubv:d 



i=3,10 



j=3,10 



c^{^I)Qri^,) + C2{^l)Q^if^) 



Cl{^,)Qf^{^i) + C2{^l)QT{^^) 



dcu I 



c^{iJi)Qr{ii)+c2{ii)Qr{ii) 



Ci(/x)Qf^(/i) + C72(;u)Qr(/^) 



ducf 



.(5) 



A basis of operators convenient for our considerations is the one in which all operators are written in 
the colour singlet form: 



Q?"'"" = {hUk){y-A){Ujdi){y-A) , Q2 



Ef^diq 

QJI = {bdi\v-A){m){v^A) , 
Qfi = -2{bq)(^s+P){Qdi)(s-P) 



{bdi)(^V~A){UjUk)(V-A) , 
{bq)(V-A){qdi)(^V-A) , 



(6) 



Q?!! = H^^i\v-A)'^q{QQ){V±A) ) 



-3eg(bq)(^s+P)iqdi)(s-P) , Qih 



2^1 



{bq){V-A){qdi)(v-A) , 
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where the subscripts {V ± A) and (S ± P) indicate the chiral structures, di = {d, s}, Ui = {u, c} and 
eq denotes the quark electric charge (e^ = 2/3, Cd = —1/3, etc.). The sum over the quarks q runs over 
the active flavours at the scale ^. Note that we use here the labeling of the operators as given in 0, ^ 
which differs from ||l|-|3| by the interchange 1^2. 

Qi and Q2 are the so-called current-current operators, Q3-6 the QCD-penguin operators and Q7-10 
the electroweak penguin operators. Ci{n) are the Wilson coefficients evaluated at ^ = 0{mb)- They 
depend generally on the renormalization scheme for the operators. For example in the HV scheme one 
has, including NLO corrections and setting rnt{mt) = 170 GeV, /x = 4.4 GeV and a^^{Mz) = 0.118 

Ci = 1.105, C2 = -0.228, C3 = 0.013, C4 = -0.029, C5 = 0.009, 
Ce = -0.033, C7/a = 0.005, Cs/a = 0.060, C9/a = -1.283, Cio/a = 0.266, ^ -* 

where a is the electromagnetic coupling constant. 

The basis in (^) differs from the one used in the NLO calculations in |Q, ^ in that some of the 
operators used there are Fierz conjugates of the ones in (^). This is the case for Q2, Qa, Qq, Qs and 
QiQ. As pointed out in |^, the Wilson coefficients evaluated in the NDR scheme (anticommuting 75 
in D ^ 4 dimensions) depend on the form of the operators and the standard Wilson coefficients in 
the NDR scheme as given in |2|, ^ cannot be used in conjunction with the basis @. On the other 
hand the HV scheme is Fierz-symmetric and the Wilson coefficients in this scheme calculated in |2|, ^ 
also apply to the basis (^). We will comment at the end of Section ^ on how the expressions for the 
effective parameters given there have to be modified when the NDR scheme with the basis of |2|, |3| is 
used. 

We observe that, unless there is some huge enhancement of their matrix elements, the contributions 
of the electroweak penguin operators and Qg are fully negligible. The operators Qg and Qiq play 
only a role in decays in which for some dynamical reasons the matrix elements of current-current 
operators and QCD-penguin operators are strongly suppressed. Generally also the contributions of 
QCD-penguin operators, especially of Q3 and Q5, are substantially smaller than those of the current- 
current operators Qi and Q2. 



3 Classification of Topologies 

Following and generalizing the discussion of ref. |^^, we classify the various contributions to the 
matrix elements of the operators Qi distinguishing the different topologies of Wick contractions 
as in figs. |l] and |2[ We have emission topologies: Disconnected Emission (DE) and Connected 
Emission (CE); annihilation topologies: Disconnected Annihilation (DA) and Connected Annihila- 
tion {CA); emission-annihilation topologies: Disconnected Emission- Annihilation (DEA) and Con- 
nected Emission- Annihilation {CEA); penguin topologies: Disconnected Penguin (DP) and Con- 
nected Penguin {CP); penguin-emission topologies: Disconnected Penguin-Emission (DPE) and Con- 
nected Penguin-Emission (CPE); penguin-annihilation topologies: Disconnected Penguin- Annihilation 
(DPA) and Connected Penguin- Annihilation (CPA); double-penguin-annihilation topologies: Discon- 
nected Double-Penguin- Annihilation (DPA) and Connected Double-Penguin- Annihilation (CPA). The 
unlabeled line corresponds to a d, n or s quark for B^, and Bs decays respectively. The dashed 
lines represent the operators. The apparently disjoint pieces in the topologies DEA, CEA, DPE, CPE, 
DPA, CPA, DPA and CPA are connected to each other by gluons or photons, which are not explicitly 
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shown. A comment about these topologies is necessary at this point. 

These special topologies in which only gluons connect the disjoint pieces are Zweig suppressed 
and are therefore naively expected to play a minor role in B decays. For this reason, they have been 
neglected in previous studies [^]. However, as we shall demonstrate in the next Section, they have 
to be included in order to define scheme and scale independent combinations of Wilson coefficients 
and matrix elements. We therefore take all of them into account in the following discussion of the 
effective scheme and scale independent parameters suitable to describe i?-decay amplitudes. We will 
later discuss some approximations that will allow us to reduce the number of parameters necessary 
to parameterize non-leptonic B decays. On the other hand, if the disjoint pieces are connected by 
photons, these contributions are automatically at least of 0{a). 

Our aim is to find a parameterization of the full contribution to the amplitude, including non- 
factorizable contributions and rescattering effects. Therefore, we will keep all the contributions defined 
above and consider them to be complex, in order to take into account final state interactions. 

The large number of different possible contributions, with different chiralities and flavour struc- 
tures, requires the introduction of fourteen flavour-dependent, scale and renormalization scheme inde- 
pendent complex parameters in order to be able to describe all decay amplitudes. The fact that the 
number of effective parameters equals the number of different topologies is accidental. This relatively 
large number of parameters is necessary if we do not want to make any specific assumption about 
non-factorizable effects and FSI. Yet, as we will see, this approach offers a transparent classification of 
various possible contributions and constitutes a good starting point for approximations which would 
reduce the number of parameters. 

In general, the topologies defined above will depend on the flavour and chiral structure, and on 
the initial and final states. In fact, if factorization were to hold, all this dependence would be taken 
into account by the factorized matrix element; however, non-factorizable contributions and FSI will 
in general have a different flavour dependence and cause for example a violation of the universality of 
the parameters af^ and introduced in ref. [^] in the framework of generalized factorization. Since 
here we take into account non-factorizable effects, we keep track of the flavour dependence. 

4 Effective Parameters — Generalities 

The matrix elements of the operators Qi{fJ.) depend in general on the renormalization scale /i and on the 
renormalization scheme for the operators. These unphysical dependences are cancelled by those present 
in the Wilson coefficients Ci{fi). Due to the mixing under renormalization this cancellation involves 
generally several operators. From the phenomenological point of view, it is desirable to identify those 
linear combinations of operator matrix elements times the corresponding Wilson coefficients which are 
both scheme and scale independent. Such combinations will define the effective parameters to be used 
in phenomenological applications. 

Before entering the details, let us make a few general comments on how the scale and scheme in- 
dependent effective parameters can be found. It turns out that the fiavour structure of the operators 
Qi and Q2, inserted in the various topologies listed above, and the known behaviour of the opera- 
tors Qi under renormalization group transformations, allow us to identify the independent effective 
parameters. 

The first scale and scheme independent combinations of Wilson coefficients and matrix elements 
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DAi{qs,q2,qi;B,Mi,M2) Qi 

CAi{q3,q2,qi;B,M,,M2) 




DEAiiqs, q2, qi] B, Mi, M2) CEA^iq^, q2, qr, B, Mi, M2) 



Figure 1: Emission, annihilation and emission-annihilation topologies of Wick contractions in the 
matrix elements of operators Qj. 
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Qi Qi ^2 




DPAiiqs, q2, qi, B, Mi, M2) CP^ife, ^2, Qi, B, Mi, M2) 



Figure 2: Penguin, penguin-emission, penguin-annihilation and double-penguin-annihilation topolo- 
gies of Wick contractions in the matrix elements of operators Qi. 
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that one can identify correspond to the emission matrix elements of the current-current operators Qi 
and Q2- Denoting by {Qi)DE and {Qi)cE the insertions of Qi into DE and CE topologies respectively, 
one finds two effective parameters 

El = Ci {Qi)de + C2 {Q2)cE , 

E2 = Ci{Qi)ce + C2{Q2)de. (8) 

We have suppressed the flavour variables for the moment. They will be given explicitly in Section |5[ 
El and E2 are generalizations of a'f^{Qi)F and a^{Q2)F in the formulation of ref. | |10| |. 

The reason why two effective parameters are needed can be found by switching off QCD effects, 
in which case Ci = 1 and C2 = 0. Dependently then on the decay channel considered, the operator 
Qi may contribute either through topology DE or through topology CE, which gives Ei and E2 
respectively. There are of course channels whose flavour structure allows both {Qi)de and {Qi)cE- 
These channels can be described by a third effective parameter which is, however, a linear combination 
of El and £'2- Consequently from the point of view of scale and scheme dependences it cannot be 
considered as a new effective parameter. 

The reason why Ei and E2 are scale and scheme independent can be understood in the following 
manner. Consider operators Qi and Q2 in the case in which all four quark flavours are different. In this 
case the penguin topologies do not contribute and also the penguin operators cannot be generated by 
QCD corrections. Choosing in addition channels in which also annihilation contributions are absent, 
we observe that Ei and E2 represent, up to CKM factors, physical amplitudes for the particular 
channels in question and as such must be scale and scheme independent. 

Next, annihilation topologies must be considered. Here we can consider channels in which both 
emission and penguin topologies as well as penguin operators do not contribute. Analogous arguments 
as given for Ei and E2 allow us to identify two new effective parameters: 

Ai = Ci{Qi)dA + C2{Q2)cA, 

A2 = Ci{Qi)ca + C2{Q2)da, (9) 

where {QijoA and {Qi)cA denote the Qj-insertions into DA and CA topologies respectively. Due to 
the flavour structure of operators Qi and Q2, Ai can only contribute to decays while A2 can only 
contribute to B^ s decays. 

It should be stressed that the Ai are independent of the Ei. Indeed, the arguments given above 
demonstrate that the cancellations of scheme and scale dependences take place separately within 
emission and annihilation topologies when only current-current operators Qi and Q2 are considered. 

The last class of non-penguin contractions that we consider corresponds to the insertion of Qi and 
Q2 into emission-annihilation topologies, denoted by DEA and CEA in fig. |^. Proceeding as above, 
we can identify two new effective parameters: 

EAi = Ci {Qi) DEA + C2 {Q2) CEA, 

EA2 = Ci{Qi)cea+C2{Q2)dea- (10) 

As in the case of Ai and A2, due to the flavour structure of Qi and Q2, EAi can only contribute to 
B^ decays while EA2 can only contribute to B^^s decays. 
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We next turn to penguin contractions of current-current operators and matrix elements of penguin 
operators. Several arguments can be given to conclude that these two types of contributions should 
be combined in order to obtain scale and scheme independent effective parameters. 

First of all let us recall that the Wilson coefficients Ci{^) and C2(/i) are independent of the 
presence of penguin operators. The same applies for the insertions of Qi and Q2 into emission and 
annihilation topologies. Consequently the scale and scheme independent effective parameters £^1, £^2) 
Ai, A2, EAi and EA2 are unaffected by the presence of penguin operators. This means that the sum 
of the remaining contributions to a physical amplitude, that is of penguin contractions of Qi and Q2 
and matrix elements of penguin operators, should be separately scheme and scale independent. 

In order to see this more clearly let us focus on 6 ^ sdd decays. Consider first TicS to be evaluated 
at a scale /Ui > rric. At this scale, we have penguin insertions of current-current operators containing 
two charm quarks {Qi^2) penguin, penguin-emission, penguin-annihilation and double-penguin- 
annihilation topologies. In addition, we have emission, annihilation and emission-annihilation matrix 
elements of penguin operators (Q3_io)- Now take instead Weflf computed at a scale ^2 < rric- the 
operators Qf^2 have disappeared, and the coefficients of operators Q^_iq have changed in such a way 
as to compensate for the absence of Qf^l- 

Now similarly to the sets (£1, £2), {Ai, A2) and {EAi, EA2) one can find four effective "penguin"- 
parameters Pi, P2, -P3 and P4. The explicit expressions for them will be derived below. Here we just 
want to relate these parameters to the penguin topologies introduced in Section |3|: 

1. Pi involves the insertions of Qi and Q2 into CP and DP topologies respectively and a particular 
set of matrix elements of QCD-penguin and electroweak penguin operators necessary for the 
cancellation of scale and scheme dependences; 

2. P2 involves the insertions of Qi and Q2 into CPE and DPE topologies respectively and a 
suitable set of matrix elements of QCD-penguin and electroweak penguin operators necessary 
for the cancellation of scale and scheme dependences; 

3. P3 involves the insertions of Qi and Q2 into CPA and DPA topologies respectively and the cor- 
responding set of matrix elements of QCD-penguin and electroweak penguin operators necessary 
for the cancellation of scale and scheme dependences; 

4. P4 involves the insertions of Qi and Q2 into CPA and DPA topologies respectively and the 
remaining matrix elements of QCD-penguin and electroweak penguin operators which have not 
been included in Pi, P2 and P3. 

In spite of the fact that in each case both disconnected and connected topologies are present only 
four effective parameters exist and not eight. This is related to the fact that the insertions of Qi into 
disconnected penguin topologies and the insertions of Q2 into connected penguin topologies vanish 
because of the fiavour structure of these operators. This should be contrasted with the insertions 
of these operators into emission and annihilation topologies, where in each case both connected and 
disconnected topologies can contribute. 

The explicit expressions for the Pi, P2, P3 and P4 parameters follows: 

5 

Pi = Ci(Qi)^cP + C2(Q2)DP + E(^2i-l((32 i-l)cE +C2i{Q2i)DE) 

4=2 
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+ J2{'^iiQi)cP + Ci{Qi)Dp)+Y.{'^^i~^iQ^^~^)cA + C2i{Q2i)DA) , (H) 
4=3 i=2 

5 

CPE + C2{Q2) DPE i-l)DE + C2i{Q2i)cE ) 

i=2 

5 10 

^{C2i-l{Q2i-l)cEA + C2i{Q2i) DEA^ + ^{Ci{Qi)cPE + Ci{Qi)DPE^ , (12) 
1=2 " ~ 

5 



+ 
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P3 — Ci{Qi)cpj^ + C2{Q2)dPA i-l)DA + C2i{Q2i)cA 

i=2 
10 

+ Y.{CiiQi)cPA + a{Qi)DPA), (13) 
i=3 

5 

P4 = Ci{Qi)^Tpj+ C2{Q2)jypj+'^(C2i-l{Q2i-l)DEA + C2i{Q2i)cEA'^ 

i=2 

(14) 



j=3 



where we have denoted by {QiY^P insertion of operator Qi in a CP topology with a c-quark 
running in the loop (this corresponds to the charming penguin of ref. ||2^), and analogously for DP, 



CPE, DPE, CPA, DPA, CPA and DPA topologies. We notice that, due to the flavour structure 
of the penguin-annihilation contributions, P3 and P4 cannot contribute to decays. Moreover P4 
contributes only to final states with two flavour neutral mesons qiqi and 9292- Similarly P2 contributes 
only to states with at least one flavour neutral meson ^2^2- 

We also stress that, in order to cancel the scheme and scale dependence of the emission contractions 
of penguin operators in the first line of the expression (12) for P2, one is forced to introduce the CPE, 
DPE, CEA and DEA topologies that also contribute to P2- Analogously, to cancel the scheme and 
scale dependence of the annihilation contractions of penguin operators in the first line of the expression 
(13) for P3, one has to consider the CPA and DPA topologies that also contribute to P3. Finally, 
when considering the insertion of penguin operators in the CEA and DEA topologies (first line of the 



expression (|lj) for P4), one has to introduce CPA and DPA Wick contractions to ensure the scheme 
and scale dependence of P4. Therefore, the Zweig-suppressed CPE, DPE, CEA, DEA, CPA, DPA, 
CPA and DPA topologies are all needed to obtain a complete scheme and scale independent result. 

The Pi, P2, P3 and P4 parameters are always accompanied by the CKM factor VtbV^^^, where 
di = d,s. Penguin-type matrix elements are also present in the part of Tieg proportional to VubV*^.. 
They correspond to penguin contractions of operators Qi,2; or, more precisely, of the differences 
^Q<kuu _ Q^iC6~j ^^2'"" — (52*'^'^) • When these combinations are inserted into penguin topologies, 
they give rise to a generalization of the GIM penguins of ref. \2G\. These scale and scheme independent 
contributions are the following: 



PfM ^ Ci[{Qi)'cp-{Qircp)+C2{{Q2)hp-{Q2)lp 

P^^^'^ = Ci(^{Qi)cPE - (Qi)c'Pb) + C2(^{Q2)dPE - {Q2)dpe 
P°^^ = Ci(^{Qi)cpA - {Qi)cPa) + C2(^{Q2)dPA - {Q2)dpA 



PfM ^ +C2((Q2)^-(Q2)^ , (15) 
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and they would vanish in the hmit of degenerate u and c. The unitarity of the CKM matrix assures 
that in a given decay Pi is always accompanied by pP™ with the same index "i". However, Pi and 
pp™ are always multiplied by different CKM factors and in order to keep the latter factors explicitly 
one has to consider separately Pi and PP™. 

It should be stressed that the eight penguin parameters introduced by us differ considerably from 
the eight penguin parameters a^^ {i = 3-10) introduced in refs. pl]]-|l3|. There a given parameter 
corresponds to a particular operator Qi and the scale and renormalization scheme dependence is 
removed by calculating the corresponding matrix elements perturbatively between external quark 
states. Such an approach suffers from gauge and infrared dependences as stressed by us in ref. |l5|] . 

These problems are absent in our approach as each of the penguin parameters Pi receives contri- 
butions from all penguin operators Q3-10 which mix under renormalization. Consequently the scale 
and scheme independence of the Pi parameters is automatically assured after the inclusion of pen- 
guin contractions of Qi,2) without the necessity for any dubious perturbative calculations of matrix 
elements. Moreover we include all possible penguin topologies, while penguin contractions in penguin- 
emission and double-penguin-annihilation topologies were not considered in the factorization approach 



in refs. |11]-|13]. The same applies to emission-annihilation topologies. 

If the operator basis is chosen so that the operators are Fierz conjugates of the ones in (^) then 
the expressions for the effective parameters given above have to be modified appropriately. Denoting 
generally by {Qi)D and {Qi)c the insertions in the disconnected and connected topologies respectively, 
the modification in question is straightforward. For each Fierz-transformed operator one simply makes 
the replacements {Qi)D ^ {Qi)c- In particular, in the case of the NDR scheme the modification 
amounts to 

{Q2i)D^ {Q2i)c, i = l-5. (16) 

Since the effective parameters are renormalization scheme independent, this transformation does not 
change their numerical values. On the other hand, in the NDR scheme the operators Q2i appear in 
the colour non-singlet form and consequently the large N counting is less transparent than in the basis 



5 Effective Parameters — Flavour Dependence 

In this Section, we discuss in detail the flavour dependence of the effective parameters defined in the 
previous Section. We give explicit expressions for the effective parameters in terms of the topologies 
defined in figures |l] and |2|. Furthermore, we provide a diagrammatic derivation of the flavour structure 
of the effective parameters. In this Section, B stands generically for the B^, B^ and Bg mesons. 

5.1 The Emission Parameters 

As we said above, the first scale- and scheme-independent combinations of Wilson coefficients and 
matrix elements that we can identify correspond to the emission matrix elements of current-current 
operators {Qi and (^2): 

Ei{qi, qj,qk; B, Mi, M2) = CiDEi{qi, qj, q^; B, Mi, M2) + C2 CE2{qi, qj,qk; B, Mi, M2) , 
E2{qi,qj,qk;B,Mi,M2) = CiCEi{qi,qj,qk] B, Mi, M2) + C2DE2{qi,qj,qk; B, Mi, M2) . (17) 



13 



In order to exhibit the flavour dependence we use here and in the fohowing the notation of figs. |l] and 
§ instead of the short-hand notation {Qi)T {T = topology) used in Section ^. If factorization held, we 
would obtain 

DEl^2{q^,qj,qk■,B,Ml,M2) = {Mi\{qjqi)^v~A)\0){M2\{bqk)(v~A)\B) , (18) 
CEi,2{quqj,qk;B,Mi,M2) = ^DE2,i{qi,qj,qk; B, Mi, M2) , (19) 

and consequently 

Ei{qi,qj,qk;B,Mi,M2) = + ^€2^ {Mi\{qjqi)(^v-A)\0){M2\(bqk)(v-A)\B) ; 

E2{quqj,qk;B,Mi,M2) = + ^Ci) (Mi|(g,gi)(v_A)|0)(M2|(6(Zfc)(y_A)|S) . (20) 

5.2 The Annihilation Parameters 

The next scale- and scheme-independent quantities corresponding to the annihilation matrix elements 
of current-current operators are given by: 

Mqi, qj,qk; B, Mi, M2) = CiDAi{qi, qj,qk; B, Ml, M2) + C2 CA2{qi,qj, qu] B, Mi, M2) , 
A2{qi,qj,qk;B,Mi,M2) = CiCAi{qi,qj,qk; B, Mi, M2) + C2DA2{qi,qj,qk; B, Mi, M2) . (21) 

If factorization held, we would obtain 

Ai{qi,qj,qk;B,Mi,M2) = (^Ci + ^^2) (MiM2|(gigfe)(y-A)|0)(0|(H(y-A)|5) ; 
A2{qi,qj,qk;B,Mi,M2) = (c2 + ^Ci^ {MiM2\iqiqk)iv-A)\0){0\(bdi)(^v~A)\B) , 
where di = d, s. Due to the flavour structure, Ai only contributes to decays and A2 to B^^s decays. 

5.3 The Emission-Annihilation Parameters 

The last two independent combinations of non-penguin topologies are given by the emission-annihilation 
matrix elements of current-current operators: 

EAi{qi, qj, qu; B, Mi, M2) = CiDEAi{qi, qj,qk; B, Mi, M2) + C2 CEA2{qi, qj,qk; B, Mi, M2) , 
EA2{qi, qj, qk] B, Mi, M2) = Ci CEAi{qi, qj,qk; B, Mi, M2) + C2DEA2{qi, qj,qk; B, Mi, M2) .(22) 

5.4 The Penguin Parameters 

Each of the parameters Pi, P2, P3 and P4 discussed in Section |^ is obtained by summing the following 
contributions, with di = {d,s}: 

a) two terms containing a given penguin contraction of Q^*'^'^ and Q'2'^'^', 



b) all the possible terms obtained by replacing the c-quark loop with the insertion of the penguin 

'3-10 



operators Q^Lio as explained below; 



c) the terms obtained by replacing Qi'2^ by Q3L10 in the penguin contractions as explained below. 
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We will now explain how the contributions b) and c) can be obtained, restricting our attention to 
QCD-penguin operators QaLg- The generalization to electroweak penguin operators is straightforward 
and will be done subsequently. A diagrammatic representation of how to obtain contribution b) is given 
in figs. P-^. Consider the insertion of a current-current operator in a penguin-like diagram in fig. |3|. 
This insertion is clearly scheme and scale dependent. As we discussed above these dependences can 
be canceled by adding the corresponding contributions of penguin operators. To this end it is useful 
to develop a diagrammatic correspondence between the penguin insertions of (5i,2 and the insertions 
of penguin operators. In fig. |3| we show the basic procedure which is applied in figs. Q-^. Wherever a 
gluon can be exchanged between the loop and a quark line, this corresponds to the diagram obtained 
by replacing the u- or c-quark loop with the insertion of penguin operators. As an example, in fig. Qwe 
show the correspondence between penguin-type contractions of and emission-type contractions 

of Qf_Q- The same is shown in fig. |5| for annihilation-type contractions of Q-jLg- In fig. ^ we show 
the correspondence between penguin-emission contractions of and emission-type contractions of 
Qs-e- The same is shown in fig. |^ for emission-annihilation contractions of Q'^'_q- In fig. |^ we show 
the relation between penguin-annihilation contractions of Qi2'^ and annihilation- type contractions of 
Q'^_Q- In fig. P we show the relation between double-penguin-annihilation contractions of Qi*!^'^ and 
annihilation-emission contractions of Q'^L^- The single gluon exchange in figs. is just a symbolic 
representation for the exchange of an arbitrary number of gluons between quark lines. 

Using the correspondence rules represented in figs. |^-|9|, we can obtain contribution b) to the scale- 
and scheme-independent combinations of penguin-like topologies. As an example, in fig. [l^ we show 
diagrammatically the construction of the b) part of the Pi contribution defined in Section ^. One starts 
by considering penguin-like contractions of current-current operators, as depicted in the first box in 
the figure, which contains the combination Ci CPi{c, dj, 92; B, Mi, M2) + C2DP2{c, di,q2] B, Mi, M2). 
Then one considers the gluon exchanges drawn in the second box, and substitutes the shaded re- 
gions with the corresponding insertions of penguin operators. The resulting third box contains the 
combination 

(CsCEsidi, q2, 92; B, Ml, M2) + C^CE^{d,,q2, ^2; B, Mi, M2) 
+CiDEi{di,q2,q2;B,Mi,M2)+CGDEG{di,q2,q2;B,Mi,M2)] 
+ (C3 CAsid^, q2,qi; B, Ml, M2) + C5 C^(d„ 92, qi;B, Mi, M2) 

+C4DA^{d,,q2,qi;B,Mi,M2) + CeDAe{d^,q2,qi;B,Mi,M2)) , (23) 



see eq. (|28|), where qi = d, u or s for B^, and Bg decays respectively. Analogously, one can build 



up part b) of P2 as shown in fig. 11. The starting point is the insertion of Qi and Q2 in penguin- 
emission diagrams, as shown in the first box in the figure. This corresponds to the combination 
Ci CPEi{c, di,q2; B, Mi, M2) -|- C2DPE2{c, di,q2; B, Mi, M2). Then, a gluon can connect the c-quark 
loop either with the qi line or with the q2 one, as shown in the second box of fig. |ll|. In the third box, 
one replaces the gluon exchange with penguin operators, following figs. ^ and 0, and one obtains the 
contribution 



C3DE3{q2, q2, df, B, Mi, M2) + C^DE^{q2, q2, d,- B, Mi, M2) 
+CiCEi{q2, q2, di- B, Mi, M2) + Ce CEe{q2,q2, d,; B, Mi,M2, 



15 




Figure 3: The basic step to relate penguin insertions of current-current operators to emission, annihi- 
lation and emission-annihilation insertions of penguin operators. Wherever a gluon can be exchanged 
between the c-quark loop in the penguin insertions and a quark line, this corresponds to the insertion 
of a penguin operator. 
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Figure 4: An example of how penguin contractions of current-current operators are related to insertions 
of penguin operators in emission topologies. First, one notices that a gluon might be exchanged 
between the c-quark loop and the q2 line in the penguin topologies (shaded region). Next, using 
the prescription in fig. |^, one replaces the gluon exchange in the shaded region by the corresponding 
insertion of a penguin operator. 
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91 





Figure 6: The same as figs. |^-|5|, for penguin-emission topologies. Tlie shaded region corresponds to a 
gluon exchange in the penguin topologies, and to the insertion of a penguin operator in the emission 
topologies. 
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Figure 7: The same as fig. ^, for another possible gluon exchange: in this case, the gluon connects the 
c-quark loop with the qi line in the penguin-emission topology. 




Figure 8: The same as fig. ^, but for penguin-annihilation topologies. 




?2 





?2 




Figure 9: The same as fig. B, but for double-penguin-annihilation topologies. 



+ (Cg CEA;{di, quq2; B, M2, Mi) + C5 CEA^{di, gi, 92; B, M2,Mi) 
+CiDEAi{di,qi,q2] B, M2, Mi) + CeDEAeidi, qi, q2; B, M2, Mi, , , 



(24) 



see eq. (p9[), where qi = d, u oi s for B^, B^ and Bg decays respectively. The same is shown in fig. 12 for 
the P3 contribution defined in Section ^. Here one starts by the insertion of current-current operators in 
penguin-annihilation diagrams (first box in the figure), which is given by Ci CPAi{c, q2, qi; B, Mi,M2)+ 
C2DPA2{c, q2,qi; B, Mi, M2). As in the previous cases, there are two possible types of gluon exchanges, 
depending on whether the gluon connects the blob to the qi or to the q2 line (second box in fig. p^). 
Substituting the gluon exchange with insertions of penguin operators, following fig. ^, we get the 
contribution in the third box, which corresponds to 

(C3DA-i{qi,q2,qi;B,Mi,M2) + C5DAr,{qi,q2,qi;B,Mi,M2) 
+C4 CA^iqi ,q2,qi;B,Mi,M2) + Ce CAe{q, ,q2,qi; B, M,, M2)) 
{C3DAs{q2,qi,q2] B, M2, Mi) + C5DAr,{q2, gi, 92; B, M2, Mi) 
+CiCA4{q2,qi,q2] B, M2, Mi) + C6C^6(g2, qi, q2; B, M2, Mi)) , (25) 



+ 



see eq. (pOD. Finally, the b) part of P4 is represented in fig. 13, The starting point here is the insertion 
of Qi and Q2 in double-penguin-annihilation diagrams, as shown in the first box in the figure. This 
corresponds to the combination CiCPAi{c,qi,q2; B, Mi, M2) + C2DPA2{c,qi,q2; B, Mi, M2). Then, a 
gluon can connect the c-quark loop either with the qi line or with the 52 one, as shown in the second 



box of fig. |l^. In the third box, one replaces the gluon exchange with penguin operators, following 
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fig. and one obtains tlie contribution 

(Ca DEA^ (g2 , (?2 , gi ; S , Ml , M2 ) + Cs DEA^ {q2,q2,qi; B, Mi, M2) 
+CiCEA^{q2, q2, qi;B, Mi, M2) + CeC^^efe, 92, qi;B, Ml, M2) 
+ (C3DEA3iqi,qi,q2; B, M2, Mi) + C5DEA5{qi,qi,q2; B, M2, Mi) 

+C4 (7E^4 (gi , 91 , 92 ; S , M2 , Ml ) + Ce C^^6 (<?! , 91 , 92 ; 5 , M2 , Ml ) ) , (26) 



see eq. (|3H). 

Concerning the c) contribution, it can be obtained by replacing the insertion of current-current 
operators in penguin topologies with the insertion of penguin operators in the same penguin topologies. 
An example is shown in fig. |l^ for the insertion of Qi'l'^ in penguin topologies, where di can be either 
a d or an s quark. The c) contribution is obtained by replacing Qi'2^ with Q'^_q, as shown in the 
figure. We remark that the Q^'_q operators contain a term with the flavour structure (bdi){didi). This 
term contributes with two different Wick contractions, depending on which di field is contracted with 
the di. This extra Wick contraction gives rise to the two diagrams in the last line of fig. 

'3-6 



The contributions of electroweak penguins to b) and c) can be obtained by replacing Qa^-e '^ith 



1=2 
5 

+ E 

1=2 



IQ using the following correspondence rules: 

Qs ^ Qg ) Q4 ^ Qw , 

Q5 ^ Q7 , Qe ^ Qs- (27) 

Applying the above considerations, we are able to identify the scale- and scheme-independent 
contributions involving penguin operators and penguin contractions: 

Pi{di,qj;B,Mi,M2) = Ci CPi(c, d„ g,; Mi, M2) + C2L'P2(c, d„ g,; 5, Mi, M2) (28) 

5 

+ ^ (C2«_ 1 CE21-1 {d, , qj ,q,;B,Mi,M2) + C21 DE21 {d, , g,- , g,- ; 5 , Mi , M2 )) 

^(C2,-iL'P2/-i(9, di, qj; B, Mi, M2) + C2iCP2iiq, di, qf,B, Mi, M2 
- q 

+C21-1 CP 21-1 {di,di,qj;B,Mi,M2) + C21 DP 21 {di,di, qy, B, Mi, M2] 

5 

+ Y.{^2l-iCA2i-i{di,qj,qk;B,Mi,M2)+C2iDA2i{di,qj,qk;B,Mi,M2)) , 
1=2 

where qk = d, u or s for B^, B^ and Bg decays respectively and q = u, d, s and c. The second and the 
last line on the r.h.s. of eq. (^8|) can be readily found by following fig. 10, The relevant quark variables 
di and qj in Pi can be identified from the penguin insertion of Qi as seen on the top of fig. The 
variable di denotes the quark line flowing into the penguin. The variable qj denotes the quark line 
attached to di on its right hand side. 

P2{d^,qj-B,Mi,M2) = CiCPEi{c,qj,di;B,Mi,M2) + C2DPE2{c,qj,d,;B,Mi,M2) (29) 

5 

+ E(C'2/-i^^2n(9i, 9j, c^*; B, Ml, M2) + C21 CE2i{qj,qj,d,; B, Mi, M2)) 
1=2 
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Figure 10: A diagrammatic representation of the correspondence between penguin-like contractions 
of current-current operators and insertions of penguin operators. See the text for details. 
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Figure 11: A diagrammatic representation of the correspondence between penguin-emission contrac- 
tions of current-current operators and insertions of penguin operators. See the text for details. 
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Figure 12: A diagrammatic representation of the correspondence between penguin-annihilation con- 
tractions of current-current operators and insertions of penguin operators. See the text for details. 
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Figure 13: A diagrammatic representation of the correspondence between double-penguin-annihilation 
contractions of current-current operators and insertions of penguin operators. See the text for details. 
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Figure 14: A diagrammatic representation of the rule to obtain contribution c). See the text for 
details. 
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+ Y.{'^2l-iCEA2i-i{di,qk,qj;B,M2,Mi)+C2iDEA2i{di,qk,qj;B,M2,Mi 



1=2 

5 



+ E 



1=2 



J2{C2i^iDPE2i^iiq, qj,di; B, Mi, M2) + C2iCPE2i{q, qj,di; B, Mi, M2 



+C21-1 CPE 21- M,qj,di-B,Mi,M2)+ C21 DPE21 {di , qj ,di; B, Mi, M2) 
where q^ = d, u oi s for B^, B^ and Bg decays respectively. The second and the third hne on the 



r.h.s. of eq. (p9|) can be found by following fig. 11. The relevant quark variables di and qj in P2 can 



be identified from the relevant penguin insertion of Qi as seen on the top of fig. 11. The variable di 
denotes the quark line flowing into the penguin. The variable qj denotes the flavour in the neutral 
"blob" representing the meson qjqj- 



P3{qj,qk; B, Ml, M2) = Ci CPAi{c, qj,qk; B, Mi, M2) + CaDP^aCc, qj,qk] B, Mi, M2) 



(30) 



^ (C2/_ 1 DA21-1 {qk,qj,qk;B,Mi,M2)+ C21 CA21 {qk , qj , qk] B , Mi, M2 
1=2 
5 

+ Yl{'^2l-iDA2i-i{qj,qk, qj; B, M2, Mi) + C21 CA2i{qj,qk, qf B, M2, Mi) 

^ (C2i^ I DPA2i^ iiq,qj,qk;B,Mi,M2)+C2i CPA21 {q, qj, qk; B , Mi, Ms 
+C21-1 CPA2i-i{di,qj,qk; B, Mi, M2) + C2iDPA2i{di, qj,qk; B, Mi, M2 



1=2 
5 

1=2 
5 

+ E 

1=2 L q 



where di = d or s for Bd and Bg decays respectively (due to the flavour structure, P3 cannot contribute 



to charged B decays). The second and third lines on the r.h.s. of eq. (|30|) can be found using fig. 12 



The relevant quark variables qj and qk in P3 can be identified from the relevant penguin insertion of 



Qi as seen on the top of fig. 12. The variable qj denotes the external quark line in the right hand side 
of the diagram. The variable qk denotes the internal quark line in this part of the diagram. 



PAiqj,qk; B, Ml, M2) = Ci CPAi{c, qj,qk; B, Mi, M2) + C2DPA2{c, qj,qk; B, Mi, M2 



(31) 



+ ^(C72z_iD^A2/-i(gi, gj, cZfc; B, Ml, M2) + C21 CEA2i{qj,qj,qk; B, Mi,M2 
1=2 
5 

+ E (C2«- 1 DEA2i^i iqk,qk,qj;B,M2,Mi) + C21 CEA21 {qk,qk,qj; B, M2 , Mi 

J2{C2l-lDPA2i-i{q, q,,qk; B, Mi, M2) + C2iCPA2i{q, q„qk; B, Mi, M2 



1=2 

5 

+ E 

1=2 



-C21-1 CPA21-1 {di ,qj,qk;B,Mi,M2) + C21 DPA21 {di, qj, qk; B , Mi, M2) 



where di = d or s for B^ and Bg decays respectively (analogously to P3, due to the flavour structure, 
P4 cannot contribute to charged B decays). The second and third lines on the r.h.s. of eq. (^) can be 



found using fig. |13[ The relevant quark variables qj and qk in P4 can be identified from the relevant 
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penguin insertion of Qi as seen on the top of fig. 13. The variable qj denotes the flavour in the upper 
neutral "blob" representing the meson qjqj. The variable qk denotes the corresponding flavour in the 
lower neutral meson (jkqk- 

If factorization held, the penguin-type matrix elements would vanish and one would be left with 
the factorized emission and annihilation matrix elements of penguin operators. 



5.5 The GIM-Penguin Parameters 

The last four effective parameters are the GIM-penguin ones, P^™, P2^^i P^™ and P^™, already 
introduced in eq. (|l5|). Their explicit flavour structure is given as follows: 



Pf™(d„(?,;S,Mi,M2) 

P^'^iqj,qk;B,Mi,M2) 
Pf™(g„(?fc;B,Mi,M2) 



+ 



+ 



+ 



= Ci[CPi{c,di,qf,B,Mi,M2) - CPi{u,d,,qf,B,Mi,M2: 
+ C2(DP2ic, di, qj;B, Mi, M2) - DP2{u, di, qj;B, Mi,M2 



Ci[CPEi{c,qj,d^;B,M^,M2) 
C2[DPE2ic,qj,df,B,Mi,M2) 
Ci[CPAi{c,qj,qk;B,Mi,M2) 
C2{DPA2ic,qj,qk;B,M^,M2) 



CTEi{u,qj,df,B,Mi,M2 
DPE2{u,qj,df,B,Mi,M2 
CPAi{u,qj,qk;B,Mi,M2: 
DPA2iu,qj,qk;B,Mi,M2) 



Ci[CPAi{c,qj,qk;B,Mi,M2) 
C2(pPA2ic,q,,qk;B,Mi,M2) 



CPAi{u,qj,qk;B,Mi,M2, 



DPA2{u,qj,qk;B,Mi,M2 



(32) 



6 Hierarchies 

In the previous Sections we have introduced the effective scheme and scale independent parameters 
suitable for the description of two-body B decays, and we have discussed in detail their flavour 
structure. In this Section, we show, using the large N approximation, that a hierarchy is expected 
between the various parameters which we previously introduced. 

Using the large N classiflcation one has in units of \/iV the following hierarchy for various topolo- 
gies: 



DE, DA 

CE, CA, DEA, CP, DPA, DPE 



CEA, DP, CPA, CPE, DPA 



CPA 



0{l/N), 
0[l/N^), 



(33) 



On the other hand the Wilson coefficients Ci_6 have in the strict 1/A^ classification the following 
dependence: 



C2 , C4 , Ce 
C3 ) C's 



0{l/N), 



(34) 
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The coefficients Cy-io are all 0{a) and consequently the large N classification of Cr-io with respect 
to Ci_6 is not useful here. In discussing the 1/N hierarchy of the effective parameters we will therefore 
not include the electroweak penguin operators. 

Two additional comments on ( ^3|) and (^) should be made. If there is an s-channel resonance, as 
in the case of annihilation contributions, there is an additional enhancement factor as in the large 
A'^ limit the width of a resonance behaves as T ~ On the other hand annihilation is suppressed 

by other dynamical factors and it is justified to ignore this additional factor of N in the following 
discussion. 

Concerning (p^), the Wilson coefficients can be enhanced by large logarithms due to the renor- 
malization group evolution from Mw down to rrif). Since In M^/m^ ~ 6, these logarithms may 
over-compensate the suppression. On the other hand as seen in (0) Ci > C2 > C4^q > 6*3^5. Thus 
except for C2 being substantially higher than C^^q the hierarchy in is roughly respected and we 
will consider it to be valid in the following discussion. 



Using (|3!^ ) and (34) we find the following large N classification of the effective parameters in units 
of VN: 



El, Ai 

E2, A2, EAi, Pi, Pf™ 
EA2, P2, P3, P^'^\ Pt"" 
Pa. Pf™ 



(35) 

0{l/N), 

0{1/N^), 

Oil/N^). 



As one can see, within the large A^ classification the annihilation amplitude Ai is expected to be larger 
than E2, which contradicts the usual expectations. Future measurements of channels dominated by 
Ai or E2 will teach us whether this hierarchy is consistent with the data. 

Next, the decay amplitudes are linear combinations of effective parameters multiplied by CKM 
factors. The hierarchy of the latter can be roughly described in terms of the Wolfenstein parameter 
A = \ Vus\ = 0.22. In particular we have: 

VusV:, = 0{X^), VtsV,l = 0{\^), V^,V:, = 0{X^), VtdV,l = 0{X^). (36) 

When calculating branching ratios one could for instance neglect VusV*^ with respect to VtsV^l 
provided this is also supported by the hierarchy in (|35|). On the other hand when studying CP- 
asymmetries it is important to keep all CKM factors and we will do so in the following. 

In tables ||-|3| we collect a large number of two-body B^, and Bs decays respectively. We list 
there the effective parameters contributing to each decay, the size of each parameter according to the 
large A^ classification and the order in A of the CKM parameters multiplying it. These tables should 
be useful in identifying the most suitable decays for the determination of the effective parameters and 
for approximations which would reduce the number of parameters. 

In the following Sections we will present explicit expressions for the amplitudes of the two-body 
B decays listed in tables In order to simplify the presentation we will omit in Sections ^ and ^: 

1. the contributions proportional to the parameters P4 and P^^™; 

2. the contributions which vanish in the SU (2) limit, unless their neglect would make a CKM factor 
disappear from the amplitude. 

The contributions omitted in Sections |^ and ^ are listed for completeness in Appendix 
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Channel 


CI. 


El 
1 


E2 
1 

N 


EA2 

1 

W 


A2 

1 

N 


Pi 

1 

N 


P2 
1 

W 


1 

w 


^GIM 
1 

N 


pGIM 
1 

w 


pGIM 
-'3 
1 

w 


^4 

1 

w 


pGIM 
1 

w 


Bd ^ D-7r+ 


A 




— 


— 


A2 


— 


— 


— 


— 


— 


— 


— 


— 


Bd ^ ^°7rO 


A 


— 


A2 


[A2] 


A2 


— 


— 


— 


— 


— 


— 


— 


— 


Bd D-K+ 


B 


A3 
























Bd D^K^ 


B 


— 


A3 






















Bd D^K^ 


B 


— 


A3 






















Bd ^ D + TT- 


B 


A3 
























Bd ^ J/^K"^ 


C 


— 


A2 


— 


— 


— 


A2 


— 


— 


A^ 


— 


— 


— 


Bd ^ 


C 


A2 


— 


— 


— 


A^ 


— 


— 


A^ 


— 


— 


— 


— 


Bd TfO^o 


D 


— 


A3 


[A3] 


A3 


A3 


[A3] 


A3 


A3 


[A3] 


A3 


[A3] 


[A3] 




D 


A3 


— 


— 


A3 


A3 


— 


A3 


A3 


— 


A3 


— 


— 




D 


A3 


— 


— 


A3 


A3 


— 


A3 


A3 


— 


A3 


— 


— 




D 


— 


A3 


A3 


— 


— 


A3 


— 


— 


A3 


— 


[A3] 


[A3] 




E 


A^ 




— 


— 


A2 




— 


A^ 




— 


— 


— 




E 




A 






\ 2 

A 


[A^] 




A 


r \ 4i 

[A*] 








Bd ^ i^V 


E 










A2 


A2 






A^ 










F 










A3 




A3 


A3 




A3 








F 






A3 






A3 






A3 




[A3] 


[A3] 


Bd^^^ 


F 














A3 






A3 


A3 


A3 


Bd ^ 


G 








A2 




















G 






A2 


A2 




















G 






A2 




















Bd ^ K+i^- 


G 








A3 






A3 






A3 






^ D+DJ 


G 








A3 






A3 






A3 








G 








A3 






A3 






A3 







Table 1: Summary of the classification of two-body Bd decays discussed in Section]^ For each channel, 
the order in the Wolfenstein parameter A of the relevant contributions is given. The contributions 
in squared brackets are of order a and they vanish in the limit of SU{2) symmetry (see Appendix 
these contributions, together with the ones in the last two columns, have been neglected in the 
analysis of Section wherever possible. 
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Channel 


Class 


El 
1 


E2 

TV 


Ai 
1 


EAi 

N 


Pi 

TV 


P2 

'W 


pGlM 

N 


dgim 
'W 




A 




A^ 














B+ D^K+ 


B 




A3 














— > U K 


r> 




\3 
A 


\3 
A 












B+ ^ L)+vrO 


B 


A3 






[A3] 












C 




A2 




A^ 




A2 




A4 




C 


A2 




A^ 




A2 




A4 






D 


A3 


A3 


[A3] 


[A3] 


[A3] 


[A3] 


[A3] 


[A3] 


B+ ^ K+ifO 


D 






A3 




A3 




A3 




5+ ^ 7r+ J/V' 


D 




A3 




A3 




A3 




A3 




D 


A3 




A3 




A3 




A3 






E 






A^ 




A^ 




A^ 






E 


A^ 


A^ 


A^ 


[A^] 


A^ 




A^ 


[A^] 


B+ K+(t) 


E 






A^ 


A^ 


A2 


A2 


A^ 


A^ 


B+ (j)-K+ 


F 








A3 




A3 




A3 


B+ ^ D+K^ 


G 






A3 












B+ ^ D+cp 


G 






A3 


A3 










B+ ^ J/V; 


G 






A3 


A3 











Table 2: Summary of the classification of two-body B^ decays discussed in Section 0. For each channel, 
the order in the Wolfenstein parameter A of the relevant contributions is given. The contributions in 
squared brackets are of order a and they vanish in the limit of SU{2) symmetry (see Appendix^); 
they have been neglected in the analysis of Section |^ wherever possible. 
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Channel 


CI. 


El 
1 


E2 

1 

N 


EA2 

1 

W 


A2 

1 

N 


Pi 

1 

N 


P2 
1 


1 


pGim 
1 


pGIM 
^2 

1 

w 


pGIM 
-'3 

1 

w 


1 

W 


pGIM 
1 




A 




- 


- 


- 


- 


- 


- 


- 


- 


- 


- 


- 


Bs ^ D^K^ 


A 


- 


A2 


- 


- 


- 


- 


- 


- 


- 


- 


- 


- 


Bs ^ D-K+ 


B 


A3 


- 


- 


A3 


- 


- 


- 


- 


- 


- 


- 


- 


Bs ^ D+K- 


B 


A3 


- 


- 


A3 


- 


- 


- 


- 


- 


- 


- 


- 


Bs ^ 


B 


- 


A3 


A3 


- 


- 


- 


- 


- 


- 


- 


- 


- 


Bs ^ Z)V 


B 


- 


A3 


A3 


- 


- 


- 


- 


- 


- 


- 


- 


- 




B 


A^ 


- 


- 


- 


- 


- 


- 


- 


- 


- 


- 


- 




B 


- 


A^ 


- 


- 


- 


- 


- 


- 


- 


- 


- 


- 


Bs ^ </.7r0 


B 


- 


A^ 


A^ 


- 


- 




- 


- 




- 




[Al 




C 


- 


A2 


A2 


- 


- 


A2 


- 


- 


A^ 


- 


A2 


A4 




C 


A2 


- 


- 


A2 


A^ 


- 


A2 


A^ 


- 


A^ 


- 


- 




D 


- 


A3 


- 


- 


A3 


[A3] 


- 


A3 


[A3] 


- 


- 


- 




D 


A3 




- 


- 


A3 




- 


A3 




- 


- 


- 


tSs J jiplX 


T~V 
U 




A 








A 






X3 
A 








Us ^ u Us 


T~V 
U 


A 








A 






X3 
A 










£>g — > XV t\ 


TT 
Cj 










A 




A 


A 




A 






Bs^(\>4> 


E 


_ 


_ 


_ 


_ 


A^ 


A^ 


A2 


A^ 


A^ 


A^ 


A2 


A^ 


Bs TT+vr^ 


E 


— 


— 


- 


A4 


- 


— 


A2 


— 


— 


A^ 


- 


- 


Bs ^ vrOvrO 


E 


- 


- 




A^ 


- 


- 


A2 


- 


- 


A^ 


[A^] 


[Al 


^ K+K- 


E 


A^ 






A^ 


A2 




A2 


A^ 




A^ 






Bs ^ i^V 


F 










A3 


A3 




A3 


A3 










G 








A2 






A2 






A^ 






^ ^0 J/V 


G 






















[A^] 


[A^] 




G 








A3 


















Bs - ^0 y 


G 






[A3] 


A3 


















Bs ^ J/V'L'O 


G 






A3 


A3 


















5, ^ J/V^L>o 


G 






A3 


A3 



















Table 3: Summary of the classification of two-body Bs decays discussed in Section^ For each channel, 
the order in the Wolfenstein parameter A of the relevant contributions is given. The contributions in 
squared brackets are of order a and they vanish in the limit of SU{2) symmetry (see Appendix |A|); 
they have been neglected in the analysis of Section ^, together with the contributions in the last two 
columns, wherever possible. 
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7 Classification of Two-Body B Decays 



Using the effective parameters defined above, we can classify the two-body B decay channels according 
to the effective parameters entering in the decay amplitude. This classification enables us to identify 
subsets of channels that, when measured, would allow us to directly extract the effective parameters 
previously defined, making no assumption about non-factorizable contributions and rescattering. We 
postpone the discussion of channels with rj and rj' in the final state to a future publication. 

As we discussed in the previous Section, we neglect here, wherever it is possible, contributions 
proportional to the parameters P4 and P^™, and contributions that vanish in the SU{2) symmetric 
limit. The expressions for the neglected terms, denoted in this Section by A^, can be found in 
Appendix 

We use the following conventions for the flavour content of mesons: 

= bu, Bd = bd, Bs = bs, vr^ = du, vr'^ = {dd — uu) , vr" = —ud, 
K+ = su, = sd, = ds, K~ = -us, D+ = sc, Dj = cs, (37) 

= dc, = cd, = —uc, = cu, cp = ss, J /il) = cc. 



which agree with refs. [18|. 



Class A decays 

Class A decays are CKM-allowed penguin-free decay channels. These are particularly interesting since 
they would allow us to extract Ei and £"2. A typical example is given by i? ^ Dtt decays. We have 
(here and in the following we give results in units of Gp/V^)'- 



A{Bd^D-^+) = VudV:f,[Ei{d,u,c;Bd,7T+,D-)+A2{c,d,u;Bd,D-,7T^ 
AiBd^D^n'') = ^^^[E2ic,u,d;Bd,D'',TT^)-A2ic,u,u;Bd,D'',Tr^))+AA{Bd^D^7 
AiB+^D\+) = VudV:k{Ei{d,u,c;B+,7r+,D'^) + E2ic,u,d;B+,D'^,7r+)) . (38) 

Other channels in class A are obtained by replacing D by D* and/or vr by p. 



Class B decays 

These are penguin-free CKM-suppressed decay channels. A typical example of class B channels is 
given by — > DK decays. We have: 

A{Bd^D-K+) = VusV:f,Ei{s,u,c;Bd,K+,D-); 
AiBd^D^K"^) = VusV:,E2{c,u,s;Bd,D'^,K^); 
A{B-^^D^K+) = VusV*i,(Ei{s,u,c;B+,K+,D^)+E2{c,u,s;B+,D'^,K+)) . (39) 

The above channels are very interesting since they are also annihilation-free and therefore would 
provide us with a cleaner measurement of emission matrix elements with respect to class A decays. 
The same holds for the analogous channels with vector mesons. 
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Other interesting class B decay channels are the following: 



AiBa^D^K^) = -VME2{u,c,s;Ba,D^,K^), 
A{B+^D''K+) = -VcsV:t,[E2{u,c,s;B+,D°,K+)+Ai{s,u,c;B+,K+,D'')), 
AiBd^Dtn-) = -VcsV:^Eiis,c,u;Bd,Dt,7r-), 



plus the corresponding ones with vector mesons. 
Claiss C decays 

Class C decays are CKM-allowed channels in which penguin contributions are present (but not dom- 
inant). Here are some typical examples: 



A{Bd ^ J/^K^) = VcsV:, E2{c, c, s; Ba, J/i>, K^) - VtsV^j, P2{s, c; B^, J/V', K^) - 

VusV:^P!i'^{s,c-Bd,J/^l^,K^); 
A{B+ ^ J/^K+) = VM Mc, c, s; B+ , J/V, K+) - VtsV:, ^2(5, c; B+, J/^, K+) - 

VusV:,[p^''-'{s,c;B+,J/i;,K+) - EAi{s,u,c; B+ , K+ , J/i^)) ; 
A{Bd ^ DtD-) = VcsV:, Ei{s, c, c; B^, Df, D") - VtsVtl P,{s, c; Ba, Df, D") - 

VusV:,P^'^{s,c;Ba,Dt,D-); 
A{B+ ^ D+^O) = VM Ei{s, c, c; B+,Dt, D") - VtsV;, Pi{s, c; B+,Dt. ^°) + 

VusV:,(A^{s,c,u;B+,Dt,D^) - P^'''{s,c;B-^,Dt,D^)) . (41) 



Class D decays 

Class D decays arc CKM-suppresscd decays in which penguin contributions are present. Well-known 
examples of this kind are i? — > tttt decays: 



V2A{Bd ^ 7r°7r°) = - {E2{n, u, d; B^, 7r°, 7r°) - A2{u, u, u; B^, 7r°, 7r°) + Pf '^(d, d; B^, 7r°, 7r°) 



A{B+ ^ I^+ttO) 




Ei{s,c,u;B+,Df,Tr^) + AAiB 



•+ 



(40) 



A{Bd ^ TT+TT-) 



+ 




)) 



+ 



Pr\d, u- Bd, 7r+, TT-) - P3«'^(«, d; B^, 7r+, Tf-)) 
VtdVtl (Pi {d, u;Bd,7T+,7r-)+ P3 (^x, d; 5^, 7r+ , vr" )) 



^(P+ ^ 7r+7r°) 





:, n, n; P+, 7r+, 7r°) E2{u, u, d; P+, 7r°, 7r+ 
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.0 



P2{d, d;B+,7r'',Tr+) - P2{d, u; S+, vr^ vr 

.0\ 



(42) 



where the quantities in square brackets in A{B~^ — > tt^tt^) vanish in the hmit of exact SU{2) symmetry. 

The presence of many different contributions with different weak phases and potentially different 
strong phases implies that it will be very difficult to extract sin 2a from the measurement of the 



asymmetry in vr+vr |26, 31 1. 

Other channels in this class are the following: 



A{Ba ^ J/V) 



+ 



E2{c, c,d; Bd,J/iJ,Tr^) 



V2 



P2{d,c; Sd, J/?/;,7r0) 



A{Bd 



D+D- 



AiB+- 
A{B+ - 
A{B+ 



vr+J/V') 



^ (P^'^^^d, c; J/V^, 7r°) - EA2{u, u, c; B^, vr", 

+ AAiBd^Tr'^J/^); 

+ VcdV:^ (Ei {d, c,c;Bd,D+,D-)+A2ic,d,c;Bd,D-,D+) 

- VudV:t,[P^'''{d, c; Bd, D+,D-) + P^'^{d, c; B^, D-,D+ 

- VtdVtl[Pi{d,c;Bd,D+,D-) + Ps{d,c;Bd,D-,D+)) ; 
+ VudV:f,[A^{d, s, u; B+,K'^, K+) - P^^d, s; B+^R"^, K 

- VtdV:^Pi{d,s-,B+,K\K+)- 

+ VM E2{c, c, d; B+, J/i;, 7r+) - VtdVtl P2{d, c; B+ , J/^,, tt+I 

- VudV:i,{P§'''{d,c;B+ ,J/ilj,^+) - EAi{d,u,c;B+ ,TT+ ,J/i,) 
+ VM El {d, c, c; 5+ , Z)+ , L)0) - VtdV:^, Pi {d, c; B+,D+,D'^) 
+ f Ai (d, c, n; 5+ , I)+ , 5° ) - Pf ™ (d, c; P+ , Z)+ , 5° 



(43) 



and the corresponding channels with vector mesons replacing the pseudoscalars. 
Class E decays 

This very interesting class consists of (charming) penguin dominated channels. These decays would 
be doubly CKM suppressed if there were no penguins; since the penguin contributions are instead 
CKM allowed, they are expected to dominate the decay amplitude. 

These modes have recently received a lot of attention after the observation at CLEO oi B ^ Kn 



decays |32]. As an example, we write down here the expression for these measured channels: 
A{Bd^K+7T-)= - VusV:,[E^{s,u,u;Bd,K+,TT-)-P^'''{s,u;Bd,K+,7T- 



A{B+ ^ K+ 



IT 



+ VMPi{s,u;Bd,K+,7r-); 

- ^^^^ (^i(s, u, u; B+, K+,7r^) + E2{u, n, s; B+,7r^,K+) 
Pf ™(s, u; B+, K+, 7r°) + Ai{s, u, u; B+,K+,tt^ 



+ 



V2 



Pi{s,u;B+,K+,TT°) + AA{B+ K^tt^); 
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A{B+^K\+)= + KsV;*b(^i(s,d,n;5+,K°,7r+)-Pf™(s,ci;5+,i^o,7r 
- Vt,V;,P,{s,d;B+,K^,7r+). 



.+ 



)) 



(44) 



Other modes in this class are the following ones: 



-^Piis, d; Ba, K^, tt") + A^(S, ^ K^tt^) ; 
VusV:j,{p^'''{s, s- Bd, cf>, K°) + ^(s, s; B^, <P, 

VtsVtl[Pi{s, s; Ba, 0, K°) + P2{s, s; B^, ; 
14.F4(pr'(s, S+, 0, K+) + Pf ^'^(5, s- B+, 0, - 




(P2(u, ^, 5; P<i, 7r°, K^) + Pf ^(s, d; P<i, K\ 7r°)) 



A{B+ K+(t)) 



Ai{s,s,u-B+,(t), K+) - EAi{s,u,s;B+,K+,(l))^ 
VtsVtl[Pi{s, s; B+, 4>, K+) + P2(s, s; B+, K+)) . 



(45) 



Class F decays 

These are CKM-suppressed pure penguin decays. Observation of these channels would also give us a 
measurement of penguin contributions. As an example, we write down the amplitudes for B K'^K^, 
B ^ (fm and B ^ (fxf) decays: 



A{Ba^K''K^)= - VudV:,{pr'\d,s-,B,,K',K^) + P^'''{d,s-,Ba,K^,K^)) 

- VM[Pi{d,s;Bd,K^,K^) + P3id,s;Ba,K^,K^)) ; 

A{Ba^cl>Tr°)= - ^^^[p^'^{d,s;Bd,cl>y) + EA2{u,u,s;Ba,7T'',cl>)) 

- y^P2{d,s;Ba,(t>,7r'') + AA{Bd^(l>7r''); 

^(P+^07r+)= - VudV:,[P^'''id,s-B+,cf>,7r+)-EAiid,u,s-B+,7r+,cl>)) 

- VtdVtlP2{d,s;B+,(l,,7T+); 

V2A{Bd ^ 00) = - VudV:hPr\s, s; Bd, 0, 0)) - VtdVMs, s; Bd, 0, 0) + AA^Bd ^ 00) .(46) 



Class G decays 

These decays proceed only via annihilations: a measurement of these channels would provide us with 
a direct determination of annihilation amplitudes. 

The best examples of class G decays are these CKM-allowed channels: 



A{Bd^DJK+) = VudV:t,A2{c,s,u;Bd,DJ,K+); 

A{Bd ^ 5° J/^) = V^V:,[A2{c, c, u; Bd, J/V, 5°) + EA2{c, u, c; Bd, D\ J/^) ; 
^(Pd^L>» = VMEA2{c,u,s;Bd,D'',^), 



(47) 
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which would allow the extraction of A2 and EA2- To be able to extract Ai one has to measure CKM 
suppressed decays such as the following: 



A{B+ - 


D+K^) 


A{B+ 




A{B+ - 


> DtJ/tl^) 



VcsV:Ms,d,c-B+,K'',D+); 

VM[Ai{s, s, c; B+, (/>, Dt) + EAi{s, c, s; B+ , 



(48) 



An example of decays proceeding through annihilations and penguin-annihilations is given by the 
following channels: 

A{Bd ^ = -VcdV^A2{c, u, c; B^, Z)°) + VMPsiu, c; B^, Z)°) 

-VudV:,{A2{u, c, u; Bd, D'>, 5°) - P^'^^iu, c; B^, 5°, D^)) ; 

A{Bd ^ DtD;) = VcdV:,A2{c,s,c;Bd,D;,Dt) - VudV:, Pt'^ic, s; B^, Df , D;) 

-VtdVtlP3{c,s;Bd,Dt,D;); 
A{Bd^K+K-) = VudV:,[-A2{u,s,u;Bd,K-,K+)+P^''^{u,s;Bd,K+,K-)) 

+VtdV,lP3{u,s;Bd,K+,K-). (49) 

The classification of Bd and B^ decays is summarized in Tables |l] and |2| respectively. 

8 Classification of Two-Body Bg Decays 

The same formalism introduced in Sections ^ and |7| can be applied to Bg decays, and a classification 
of the various decay channels can be made according to the CKM structure and to the effective 
parameters entering the amplitude. 

As in the previous Section, we present here the results obtained by neglecting terms that vanish in 
the limit of exact SU{2) symmetry and discarding the effective parameters P4 and ^4^'^, wherever it 
is justified to do so. The neglected terms, denoted in this Section by A^, can be found in Appendix 

i 

Class A decays 

Class A decays are CKM-allowed penguin-free decay channels. We have: 



AiBg^D;7T+) = VudV:,E,id,u,c;Bs,7T+,D;); 

A{Bs^D^K^) = VudV:,E2{c,u,d;Bs,D^,K''). (50) 



Other channels in class A are obtained by replacing pseudoscalar mesons with vector ones. 
Class B decays 

These are penguin-free CKM-suppressed decay channels. We have: 

A{Bs^D;K+) = VusV:,[e,{s,u,c;Bs,K+,D;) + A2{c,s,u;Bs,D;,K+)), 
A{Bs^DtK~) = -VcsV:,[Eiis,c,u;Bs,Dt,K-) + A2{u,s,c;Bs,K-,Dt)) , 
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A{Bs^D"^) = -VcsV:,,(E2{u,c,s;Bs,D^,^) + EA2{u,c,s;Bs,D'',^)). (51) 

Other class B decays are the following doubly-CKM-suppressed transitions: 

A{B,^D+K-) = -VMEi{d,c,u;B„D+,K-); 
A{Bs^D^K^) = -V,dV:,,E2{u,c,d;Bs,D'>,K^); 

^(S,^7rV)= - ^^(£;2(^x,^x,s;5„7^^0) + ^A2(u,u,s;S„7^^</.)) 

^*-^t6[p2(s,d;-B„7r°,(/.)-P2(s,u;5„7r0,(/.)l + A^(B, ^ ttV) , (52) 



V2 

plus the corresponding ones with vector mesons. The term in square brackets in A{Bs — > T^^(p) vanishes 
in the limit of exact isospin symmetry. 

Class C decays 

Class C decays arc CKM-allowcd channels in which penguin contributions are present (but not dom- 
inant). Here are some typical examples: 

A{B, ^ DtDj) = VM[e,{s, c, c; B^, Df, DJ) + ^2(0, s, c; S„ DJ,Dt)) - 

VtsV;^, (Pi {s, c; Bs , Dt ,D-) + P3 (c, B,, Z)+ , I); )) - 
VusV:^{P?'''{s, c- Bs, Dt, D;) + P3°™(c, s; B,,Dt, D;)) ; 

A{Bs ^ J/i^cP) = VcsV:,,[E2{c, c, s; B^, J/^j, 0) + ^^2(0, c, s; S„ J/V', </.)) 

- VtsV;, P2{s, c; B,, J/^, 0) - VusV:^ P^'^'is, c; S„ J/V, 0) 

+ AA{Bs^ J/ip(l)). (53) 

Class D decays 

Class D decays are CKM-suppressed decays in which penguin contributions are present: 

A{B,^Tr^K^)= - y^(E2{u,u,d;Bsy,K^)+P^'''{d,d;Bsy,K^)) 

- ^?^Pi(d,d;i?„7r°,i?°) + A^(i?3-^7r°i^0); 

AiB,^7T+K-)= - V^dV:i,{E^{d,u,u;B„Tr+,K-)-P^'''{d,u-Bs,Tr+,K-)) 
+ VtdV*t,Pi{d,u;Bs,TT+,K-); 
A{Bs ^ D+DJ) = VM Md, c, c; Bs,D+, DJ) - VtaV^, Pi{d, c; Bs,D+, DJ) - 

V^dV:t,P^'''{d,c;B,,D+,DJ); 
A{Bs ^ J/V) = VM E2{c, c, d; Bs,J/^, K^) - VuV*, P2{d, c; B„ J/^, K^) - 

VudV:, P^'^^id, c; Bs, J/V, ^°) . (54) 
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Class E decays 

This class contains channels in which penguin contributions are CKM-allowed, while non-penguin ones 

are either absent: 

- VtsVtl[Pi{s, d; Bs,K'', ^0) + P^{d, s- Bs,K\ K^)) ; 

V2A{B, ^(t>(t>)= - K,F4(2Pf '^(s, s- B„4>, cf>) + 2P2°i^(s, s; B„ ct>, cf>) + P^'^'is, s; B„(l>, 0)) 

- VtsVtl{2Pi{s, s; Bs, cP, </>) + 2P2{s, s; B^, </)) + Psis, s; B^, <P, <P)) 

+ AA{Bs ^ <P<P) ; (55) 

or doubly CKM-suppressed: 

A{Bs ^ TT+TT-) = + K,y4(P3°™(«,d;B„7r+,7r-)-A2(u,d,u;5„7r-,7r+)) 
+ VtsVtlP3{u,d;Bs,TT+,7r-); 
V2^(P,^7rV)= - VusV:,{^P,^'^{u,u;Bsyy) + lpt^{d,d;B,yy) 

-A2{u,u,u;B,,Tr^,TT^)^ 
- Vt,V,*,(^P3(«,t^;S„7rO,7rO) + ^P3(d, d; B„ tt^, Tr^)) + A^(P, ^ ttV) ; 
A{B,^K+K-)= - VusV:t,[Ei{s,u,u;Bs,K+,K-)+A2iu,s,u;Bs,K-,K+)- 

Pf u; Bs,K+, K-) - P^'^'iu, s; Bs,K+, R-)) 
+ VtsV^t,(Pi{s,u;Bs,K+,K-) + Ps{u,s;Bs,K+,K-)). (56) 

Class F decays 

The following channel is an example of a CKM-suppressed, pure penguin decay: 
A{Bs ^ ^V) = - VtdVtl(Pi{d, s; Bs,K^, 0) + P2{d, s; B^,^, K^)) 

- KdV;*,(pf '-(d, s; Bs,K\ 0) + Pf -(d, s- BsA. ^°)) ■ (57) 

Class G decays 

The following channels are CKM-allowed decays that proceed only through annihilation and penguin- 
annihilation: 

A{Bs ^ D-D+) = VcsV*,, A2{c, d, c; Bs,D~,D+) - K.^^ P^'^^id, c; Bs,D-,D+) 

-VtsVtlP3id,c;B„D-,D+); 
A{Bs ^ D^'D^) = -VMMc, u, c; P„ I?") + VtsV^lPsiu, c; P„ D^, Z)0) 

+K,F„*;,(n°™(^,c;P„P>o,L>°) -^2(u,c,u;P„L>o,L>°)) . (58) 

Bg — 7r°J/V' is an annihilation channel dominated by electroweak (isospin-breaking) effects: 
^(P, ^ 7r° J/^) = + ^^^^\EA2{c,c,d;Bs,J/i^,n'^)- EA2ic,c,u;Bs,J/i;,Tr^) 
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- ^^^EA2iu,u,c;Bs,7r'>,J/i;) + AA{Bs^Tr'^J/ij), (59) 
v2 

where the quantities in square brackets vanish in the SU{2) limit. 

Here are some examples of CKM-suppressed pure annihilation decays: 

: -V^sVuhA2{u,d,c;B,,T:- ,D+); 

: VMA2{c.d,u-B,,D-,n+)- 

■- -VcsV:f,[A2{u,c,c-Bs,D^,J/i;) + EA2{u,c,c;Bs,D'',J/i;)) ; 

: VusV:,[A2{c,c,u;Bs,J/^P,D'') + EA2{c,u,c;Bs,D^,J/i,)) ; 

: Y^A2{u, u, c; Bs,7r^,D^) + A^(5, ^ tt^D^) ; 

: _]:if£^^2(c,u,u;5s,5°,7r°) + A^(S, ^7r°L'°). (60) 
V2 

The classification of Bg decays is summarized in Table ^. 

9 Strategies for the experimental determination of the effective pa- 
rameters 

In this Section we would like to make general comments on the determination of the effective pa- 
rameters from the data. A detailed numerical analysis will be presented in a subsequent work. The 
inspection of tables |^-^ allows to identify the most suitable decays for this determination. We have: 

i) \Ei\ can be best determined from B^ D^K^, Df-K^, B^ Dfir^, Bg D~7r~^; 

ii) 1^21 can be extracted from Bd D^K^, D^K^ and Bg D^K^; 

iii) I All can be determined from B^ D^K^; 

iv) IA2I can be taken from Bg tt- D+ , ir+D-, tt^D^, tt^D'^; 

v) \EA2\ can be extracted from B^ — > D^(p; 

vi) |Pi| can be determined from B^ K^n" , K^tt^ and B^ K^tt^, K^n^; 

vii) IP3I can be extracted from B^ cfxj) and Bg tt^tt". 

In all these cases the parameter in question dominate a given decay or constitute the only contribution. 
An independent determination of EAi and P2 is harder: 

viii) EAi could play some role in B~^ — > vr"*' J/V', (p'^'^i E>f(j), DfJ/'ip; 

ix) P2 could be non-negligible in B~^ — > J/iJjK~^, (j)TT~^ and in particular in Bd — > J/t/jK^, (jyir^, 
Bg — > J/ipcj) and (fxf). On the other hand, this contribution being 0{1/N'^), it appears that this 
determination will only be possible by comparing various decays. 
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The determination of pP^^ is complicated by the fact that these contributions are often suppressed 
by 0{X'^) relatively to Pi contributions. However, in situations in which pP™ are multiplied by CKM 
factors with large complex phases, as VusV*^^, CP-asymmetries could be useful in this respect. 

Finally we do not think that the parameters P4 and P^™ will be determined in the near future, 
as they are expected to be very small and one would need very precise data to investigate their effect. 

When data will be available on most of the channels discussed in Sections ^ and ^, it will be 
possible to study the flavour dependence of the effective parameters and to learn more on final state 
interactions by extracting the phases of the effective parameters from the data. 



10 Comparison with the Diagrammatic Approach 



We will now compare our approach with the diagrammatic approach of refs. |17, in particular 
with the formulation given in [l^ ]. 

The main difference between these two approaches is the following one. Whereas the diagrammatic 
approach of [|^] is formulated in terms of diagrams with full W, Z and top quark exchanges, the 
approach presented here is formulated directly in terms of local operators, the basic objects of the 
effective theory. The main advantages of formulating non-leptonic decays in terms of operators are as 
follows: 

i) The effective parameters introduced in Section ^ are directly expressed in terms of matrix el- 
ements of local operators. Therefore they can be in principle calculated in QCD by means of 
suitable non-perturbative methods. Consequently the comparison of phenomenologically ex- 
tracted effective parameters from the data with the values calculated in QCD may offer some 
useful tests and teach us something about strong interactions. Such tests are clearly not possible 



in the approach of |18] as no prescription is given on how the phenomenological parameters in 



this approach could be calculated in QCD. 

ii) The inclusion of QCD perturbative corrections can be consistently performed by using the NLO 
Wilson coefficients in HeS- 

iii) The issues of renormalization scheme dependences, non-factorizable contributions, flavour de- 
pendences and in particular the important issue of final state interactions can be addressed 
transparently in the operator approach, whereas the diagrammatic approach can be in this 
context sometimes even misleading. 

Several of these advantages cannot be fully appreciated yet in view of the limitations of present 
non-perturbative methods, but this may improve in the future. 



In spite of these basic differences between our approach and the one in refs. |17, 18 1, it is possible 
to establish some connections between them and to show explicitly where they differ from each other. 

In the approach of ref. |]l8[ the basic parameters for strangeness-preserving decays are the "Tree" 
(colour favoured) amplitude T, the "colour suppressed" amplitude C, the "penguin" amplitude P, the 
"exchange" amplitude E, the "annihilation" amplitude A and the "penguin annihilation" amplitude 
PA. For strangeness-changing decays one has T' , C, P' , E' , A', PA' . If Z-penguins are taken into 
account one introduces in addition the "colour-allowed" Z-penguin P^w and the "colour-suppressed" 
Z-penguin P^^. Similarly P^^ and P^ are introduced for strangeness-changing decays. All these am- 
plitudes include the relevant CKM factors, whereas in our approach the CKM factors are not included 
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in the effective parameters. While this distinction is important for phenomenological applications, we 
will omit the CKM parameters in the discussion below, as far as possible. 



Let us investigate how the parameters of ref. |18| are related to the parameters introduced in 
Section ^ It is sufficient to consider "unprimed" amplitudes. The discussion of "primed" amplitudes 
is completely analogous. 

As stated above, the amplitudes T, C, P, E, A, PA, Pew and P^ are defined in terms of 
diagrams in the full theory which contain explicit W, Z and top propagators, whereas our approach 
is formulated in terms of operators and diagrams in the effective theory. The connection between 



these two approaches can be established by noting that in all diagrams of ref. [18| a PV^ propagator is 
present, whereas in our approach the operator Qi contributes to all the effective parameters. Since a 
1/F-exchange between two quark lines is represented in the effective theory by Qi this is not surprising. 

Now, the situation is of course more complicated as QCD corrections to a VF-exchange generate 
other operators for which no diagrams exist in the diagrammatic approach of ref. In the limit 

as = 0, a = we have, however, Ci = 1 and Cj = (i 7^ 1). Consequently in this particular limit 
the effective parameters in our approach are entirely given in terms of matrix elements of Qi and the 
correspondence between our approach and the one of ref. [^] is easier to establish. One can then think 
that when QCD and QED corrections are included, the contributions of the operators Qi {i 7^ 1) are 
added properly to those of Qi so that the effective parameters are scale and renormalization scheme 
independent. 

Proceeding in this manner it is easy to establish first the following correspondence 

T ^ El, C ^ E2, A^ Ai, E ^ A2. (61) 

The Zweig-suppressed contributions represented in our approach by EAi and EA2 have not been 
taken into account in [18|, although it is straightforward to draw the corresponding diagrams. In the 



case of charmless i?-decays to two pseudoscalars, considered in [^], it is very plausible that these 
contributions can be neglected as seen in tables |l[-|3[ On the other hand EAi could be important for 
— > (j)TT~^ and in particular in decays with charm in the final state such as tt+J/^/;, Dfcj) a-iid 

DfJ/tp. EA2 is fully responsible for D^cj) and could be significant in B^ — > -k^J/iI), Bs — > D^4>, 

The case of P, PA, Pew and P|^ is more involved. If one sets a = we have roughly speaking 

P ^ Pi, PA^Ps, (62) 

but this correspondence is a bit oversimplified and requires some explanation. It is sufficient to discuss 
only the first relation. 
Let us write P as 

P = VudV:, {Pu - Pc) + VtdVtl {Pt - Pc) , (63) 



where we have used the unitarity of the CKM matrix. In the language of ref. [1^ P„ , Pc and Pt denote 
QCD penguin diagrams with internal u, c and t exchanges. As discussed already in p^ , |27| , in the 
operator approach Pt is represented by the contributions of the QCD penguin operators Q3-6 whereas 
Pu and Pc by the matrix elements {Qi)cp i{Q2)^p) and {Qi)cp HQ'^Ydp) respectively. Thus we can 
establish the relation 

Pc-Pt^Pi, Pc-Pu^ Pi'"^ . (64) 
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At this point, following [gj], it should be emphasized that whereas Pi and P^™ are /i and renormal- 
ization scheme independent, this is not the case for P„, Pc and Pt- Consequently while P^™ could 
possibly be neglected with respect to Pi, the neglect of Pc with respect to Pt or of Pu with respect to 
Pc would automatically introduce unphysical scheme dependences. 

In this context we would like to recall that the impact of Pc and P„ on the extraction of CKM- 
phases has been investigated for the first time in In the operator language, the importance of the 
charm contribution to Pi ("charming penguins") in connection with the CLEO data on P ^ Ktt has 
been stressed in pq] and analyzed in detail in p3] . Finally the role of Pu in connection with final state 



interactions in B decays has been pointed out in |2J] and analyzed subsequently in [20|-|2^], |25]. 

Let us now discuss the issue of electroweak penguin contributions which in the approach of ref. 
are represented by Pew and P|^. It is sufficient to discuss Pew only. As in the case of P one can 
write 

Pew = VudVub i^EW - Pew) + ^id^tfe (^ew " ^ew) i (65) 



where in the language of ref. |24, 27, ^ Pewj -^ew -^ew denote Z-penguin diagrams with internal 
u, c and t exchanges. Strictly speaking P^^ have to include also the box contributions to make them 
gauge independent. 

In our approach the electroweak penguin contributions are represented by the electroweak penguin 
operators Qt-iO) which are included in Pi and in the other penguin parameters Pj. The point is that 
the contributions of Qj-io are by themselves scale and renormalization scheme dependent and have 
to be considered simultaneously with other operators to obtain physical scheme independent results. 
Consequently when electroweak penguin contributions and generally 0{a) effects are included the 
correspondences (^) generalize to 

Pc-Pt+ Pe'w - PL ^Pl, Pc-Pu + Pew " ^ew P?'"" , (66) 

with 

Pc + P^^^Ci{QiYcp + C2{Q2)hp, Pu + P^^^Ci{Qircp + C2{Q2)lp. (67) 

Here {QiYcp includes the insertion of Qi in the CP penguin topology both with gluon exchanges and 
with a single photon exchange. It can be considered as a sum of QCD charming penguins and QED 
charming penguins. Similar comments apply to {Q2)f)p, (Qi)cp ^^'^ {Q2)^p- The QED charming 
penguins as well as the QED u-penguins have been identified in |27| but neglected with respect to top 



penguins represented in our approach dominantly by the contributions of the operators Qg and Qio- 
Note that Ci^2 = C'(l) whereas (Qi)c'p {Q2)^dp with a single photon exchange are 0{a). On the 
other hand Cg and Cio are 0{a) but {Qg)cE and {Qio)de contributing to Pi in ([Tl[) are 0{1). Thus 
from the point of view of an expansion in a the QED c- and u-penguin insertions of Qi and Q2 are 
of the same order as the Qg and the Qio contributions. In fact, as we stated above, they have to be 
both included in order to obtain scheme independent results. On the other hand one could argue that 
the strong enhancement of Cg^io through the large top quark mass and the suppression of (Qi)c'p 
and (Q2)£)p through l/lGvr^ factors present in a perturbative evaluation of these matrix elements 
makes the neglect of QED c- and u-penguin insertions of Qi and Q2 plausible. This is supported to 
some extent by the very weak scheme dependence of Cg and Ciq. Still one should keep in mind that 
perturbative arguments for the smallness of 0{a) corrections to {QiYcp {Q2)'dp apply 
and the relevance of these contributions has been possibly underestimated in the literature so far. 
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The Zweig-suppressed penguin contributions represented in our approach by P2 and P4 have not 
been taken into account in As in the case of EAi and EA2 their role in charmless B decays to 
two pseudoscalars is expected to be very small. As seen in tables Q-^, P2 could play some role in a 
number of decays with charm in the final state, in (j)Tr^, K^cj) ^'^^ analogous channels in B^ and 

Bs decays. P4 being doubly Zweig-suppressed is most probably negligible in all decays. 



Finally we would like to comment on the parameters P, T, -P^, C etc. used in refs. |24, 27, 
These parameters should not be confused with the ones discussed above. They have been introduced 
in order to make the analysis of the extraction of the angle 7 from B — > Ktt decays more transparent. 
As an example we show how the parameters P, T and defined through 



A{B+ TT+K") = P , AiBd TT-K+) = -(P + T + Pj 
are given in terms of the effective parameters introduced in Sections |^ and ^. One has 

P = VusV:,[Ai{s,d,u;B+,K'^,TT+)-P^'''{s,d-B+,K^,TT+))-VtsVtlPi{s,d-B+,K^,TT+), 



(68) 



T 



VusV*A Ei{s, u, u; Bd, A'+, vr" 



+ 



pGIM| 



s,d]B+,K^,Ti^ 



Pi{s,u;Bd,K+,TT ; 
- Ai{s,d,u;B+,K^,7T+) 
-Pf™(s,n;P,,K+,7r-) 



Pi(s, d; B+, i^o, 7r+) - Pi(s, n; B^, vr" 



(69) 



where A and A are the Wolfenstein parameters. 

Analogous expressions can be found for other B — > Kn decays. 



11 Summary 

In the present paper we have proposed a general framework for analyzing non-leptonic two-body 
P-decays which combines the operator language with the diagrammatic language. Following and 
generalizing the discussion of ref. [^ ] we have classified the contributions to the matrix elements of 
the relevant operators in terms of different topologies of Wick contractions. Subsequently we have 
introduced a set of effective parameters which are both renormalization scale and renormalization 
scheme independent. As such they are convenient for phenomenological applications. 

On the other hand, being linear combinations of Wilson coefficients and particular Wick contrac- 
tions of local operators, these effective parameters are in principle calculable in QCD. This feature 



distinguishes our approach from the diagrammatic approach of refs. |17| , 18 1 in which no reference 
to local operators is made and no prescription for the calculation of the corresponding parameters is 
given. 

The formulation given here allows to describe in general terms the flavour dependence of non- 
leptonic two-body decays including non-factorizable contributions and final state interactions. It is 
therefore particularly useful for a general model-independent study of CP violation in B decays. 

In the present paper we did not use any symmetry arguments like SU{2) or SU (3) flavour symme- 
tries. On the other hand we have used the expansion to indicate a possible hierarchy among the 
effective parameters. In particular we have included in our discussion a number of Zweig-suppressed 
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topologies, which have not been discussed in the hterature. We have shown that these topologies have 
to be included in order to obtain a consistent description of non-leptonic decays with respect to scale 
and renormalization scheme dependences. While such topologies are suppressed in the large N limit, 
the role of Zweig-suppressed contributions in non-leptonic decays is an interesting and important issue, 
which requires further theoretical and phenomenological investigations. 

As a preparation for phenomenological applications of our formalism we have presented a clas- 
sification of two-body B decay channels according to the effective parameters entering in the decay 
amplitudes. This classification enabled us to identify subsets of channels that, when measured, would 
allow a direct determination of the effective parameters making no assumption about non-factorizable 
contributions and rescattering. The execution of this program is deferred to a subsequent work. 

The approach presented here allows for a phenomenological description of non-leptonic decays, 
once a large number of channels have been measured. Such a description may teach us about the 
role of non-factorizable contributions and about the flavour structure of non-leptonic B decays. One 
should hope that in this manner some regularities will be found. However, without a dynamical input 
the formulation presented so far can be considered only as the most suitable language to describe 
non-leptonic decays in a manner consistent with QCD. In order to be predictive some additional input 
involving symmetry arguments and dynamical assumptions is needed. We will return to these issues 
in a subsequent publication. 



Acknowledgments 

We would like to thank Stefan Bosch for invaluable comments on the manuscript and for checking 
many equations. We also thank A. Khodjamirian for discussions. 



A Remaining Contributions 

We collect here the contributions that were neglected in the analyses of Sections and ^. These 
include the terms proportional to P4 and ^4^'*^, and some of the terms vanishing in the limit of 
SU{2). The SU{2) breaking effects come on the one hand from the matrix elements of electroweak 
penguin operators, and on the other hand from 0(a) effects in the matrix elements of operators 
Qi-6- We stress that, for consistency, if one takes into account the effects of electroweak penguin 
operators, one should also include 0{a) effects in the matrix elements of operators Qi-q. Therefore, 
in this case operators Q3-6 can also contribute to AI = 1 (or 3/2) transitions, due to 0{a) effects 
in the matrix elements. These contributions will be, for example, of the form Pi{d,u; B~^ ,7r~^ ,Tr^) — 
Pi{d,d;B'^,TT^,ir^), vanishing in the SU{2) symmetric limit, and will be proportional to a. Thus, 
while the individual Pi penguin parameters are of order a^, differences of penguin parameters vanishing 
in the SU{2) symmetric limit are at least of order a. 

Class A decays 

We have (here and in the following we give results in units of Gp/^/2): 

AA{Bd^D\'>) = ^^^[EA2{c,u,d-Bd,D'>,7r^)- EA2{c,u,u;Bd,D'>,Tr'>)]. (70) 
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Class B decays 

We have: 



AA{B+ ^ L>+7r°) 



V2 
V2 



EAi{s, c, u; B+,D^, 7r°) - EAi{s, c, d; B+,D+, 7r°) 



+ 



V2 



P4(ci, s; Bs, 7r°, ^) - P^iu, s; Bs,n°, < 



Class C decays 

We have: 



AAiBs ^ J/iPct)) = -VtsVtlPi (c, s-Bs,J/^,4>)- VCbPr™ (c, s; P„ J/V', <t>) . 



Claiss D decays 

We have: 



0„0\ 



+ 



(^[^^2(^^, d; Bd, TfO, 7r°) - E^2(u, u, u; 5^, 7r°, vr") 
'^Pf-(^,^x;Pd,7rO,7rO) - Pf d; S^, tt^, tt^) + lp4°™(d, d; P^, 



Pf -(d, u; Bd, ttO, 7r°) - Pf -(d, d; P^, 7^^ tt' 



^0 



A^(Prf ^ TfO J/V-) = + 



VidF/b ( n; Pd, ttO, ttO) - P4(ii, d; P^, tt^, ttO) + ^P4(d, d; P^, 



P2(d, n; Pd, 7r°, Tr") - P2{d, d; B^, 7r°, tt^) 



+ 



VudVub 

VtdVl 



A^(P+ ^ 7r+7r°) = + 



V2 
V2 



P^2(c, c, ci; Prf, J/V', 7r°) - EA2{c, c, u; Ba, J/^p, vr" 
P4°™(d, c; Pd, 7r°, J/iP) - P4°™(«, c; P^, 7r°, J/V) 
P4(^x, c; Pd, 7r°, J/V') - P4(d, c; P^, 7r°, J/V') 
yli(d,d,u;P+,7r^7r+) - u; P+, 7r+, tt" 



+ 



Pyli(d, n, cZ; P+,7r+,7r°) - P^i(cZ, -u, -u; P+, 7r+, 7r°) 
Pf '^(d, d; P+, 7r°, 7r+) - Pf n; P+, 7r+, tt")] + 

P2°'"(d, u; B+y, 7r+) - P2°'-(d, d; B+y, 7r+) 
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V2 



Pa {d, d; P„ 7r°, - P2 (d, n; P., 7r°, 7?°) 



(73) 



Class E decays 

Wc have: 
AA{Ba^K\^)= - 

AA{B+ ^ K+TT^) = - 



VtsVl r 



V2 
V2 



P,^'^{s, d- B,, ttO, i^O) - P,^'^{s, u; B^, 7r',K' 
P2{s, d; Ba, 7r°, if^) - ^2(5, u; B^, 7r°, if °) 

EAi {s, u, u; B+, K+,tt°) - EAi (s, u, d; B+,K+, 7r°) 

P^'^^is, d- B+y, K+) - P^'^is, u; B+, 7r°, K+) 

P2is,d;B+,Tr^,K+)-P2{s,u;B+,n^,K^ 



+ 



AA{Bs - 



V2 



us ^ub 



-P4 (■u,M;i^s,7r ,7r ) + -P4 [d,d; Bs,n ,Tr ) - 



P4°'^'(u,(i;P„7r",7r") + E^2(«, 5., vr", tt") - P^2(^^, 5., 7r°, 7r°^ 



^P4(«, u; P„ 7r°, 7r°) + ^P4(d, d; P„ 7r°, 7r°) - P4(«, d; P„ 7r°, 7r°) 



(74) 



Class F decays 

We have: 

AA{Bd ^ <^7rO) 



= + 



+ 



VtdV. 



tb 



V2 
V2 



Pi{u, s; Bd, 7r°, (p) - Pi{d, s; Ba, 7r°, ( 



P4°'-(u, s; Bd, 7r°, 0) - P4°'-(d, 5; Ba, tt", 0) 



A^(5<i^#)= - VuV,lP^is,s;Ba,(l>,(l>)-V^dV:,Pt''is,s;Ba,(t^,4>). 



(75) 



Class G decays 

We have: 

AA{Bs ^ 7r° J/V-) = 

AAiB, ^ 7r°L>°) 

A^(P, ^ 7rOL>0) 



V2 



'P^'id, c; P„ 7r°, J/V) - P4°™(n, c; P„ 7r°, J/V) 
E^2 («, c, u; 5„ , 7r°) - EA2 {u, c, d; P, , , 7r°) 
EA2{c,u,u;Bs,D'^,7r^) - P^2(c,tt,(i; Ps,P'°,7r° 



(76) 
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